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Abstract 

Fomalont and Kopeikin have recently succeeded in measuring the 
velocity-dependent component of the Shapiro time delay of light from 
a quasar passing behind Jupiter. While there is general agreement that 
this observation tests the gravitomagnetic properties of the gravitational 
field, a controversy has emerged over the question of whether the results 
depend on the speed of light, c, or the speed of gravity, Cg. By analyzing 
the Shapiro time delay in a set of "preferred frame" models, I demon- 
strate that this question is ill-posed: the distinction can only be made 
in the context of a class of theories in which c 7^ c^, and the answer 
then depends on the specific class of theories one chooses. It remains 
true, however, that for a large class of theories "close enough" to general 
relativity, the leading contribution to the time delay depends on c and 
not Cg] within this class, observations are thus not yet accurate enough 
to measure the speed of gravity. 



* email: carlip@physics.ucdavis.edu 



1 Introduction 



In a remarkable experiment, Fomalont and Kopeikin [1] recently observed the Shapiro time 
delay of hght from the quasar QSO J0842+1835 as its image passed within 3.7 arcminutes 
of Jupiter. Using VLBI, they succeeded in measuring not only the first-order effects, but the 
dependence on Jupiter's velocity as well. The results agree well with the predictions of general 
relativity, and there is general agreement that the observation has tested the gravitomagnetic 
effect of Jupiter's motion. Within the framework of general relativity, this is essentially 
equivalent to a confirmation of the Lorentz transformation properties of the gravitational 
field. A controversy has emerged, however, over the question of whether the results depend 
on the speed of gravity, Cg — in which case the observation can be construed as a measurement 
of Cg as well — or only on the speed of light, c [2-9]. 

More precisely, the delay predicted by standard general relativity due to a gravitating 
body of mass rua and velocity Va takes the form 



where is the separation vector between the gravitating body and the observer and k is the 
unit vector in the direction of the incoming light ray. The question at issue is whether c in 
eqn. (1.1) should be understood as c or Cg, or some combination — perhaps even a different 
combination in each of its appearances. 

Let us begin with two simple observations. First, within standard general relativity, c and 
Cg are identical: gravity propagates along null geodesies [10], and c and Cg are determined 
by exactly the same metric information. This equality is deeply embedded in the Einstein 
field equations, since, as is well known, the field equations themselves determine the motion 
of the sources [11]. One thus cannot simply declare c and Cg to be two separate parameters 
in general relativity; to distinguish the two speeds, one must instead look at a broader class 
of "comparison theories" in which c to determine the separate dependencies of the time 
delay on the two speeds. Second, it is not obvious that the answer will be independent of the 
class of "comparison theories" one chooses. Indeed, Will's PPN formalism [3] and Kopeikin's 
alternative PPN formalism [4, 5] may be seen in part as two different choices that yield two 
different results. 

In this paper, I evaluate the Shapiro time delay in a class of preferred frame theories 
described by Jacobson and Mattingly [12] (sec also [13]). While these are by no means the most 
general models with c 7^ c^, they are sufficiently general to clearly show that the parameter 
c in (1.1) is model-dependent, and can depend in complicated ways on c and Cg. The debate 
over whether Fomalont and Kopeikin have measured the speed of gravity is thus misguided: 
the interpretation of their measurement depends on a choice of "comparison theory," and 
cannot be narrowed down unless that choice is strongly constrained by other means. 




with K = k — -k X (va X k). 
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I also show, however, that for theories in the class I consider — and probably for a consider- 
ably more general class of theories that are sufficiently "close" to general relativity — the factor 
of c that appears inside the logarithm of (1.1) is the speed of light. Since for ordinary general 
relativity this term gives the dominant contribution to the Fomalont-Kopeikin observation, it 
is arguably true that observations are not yet sensitive enough to determine Cg. 

2 A Class of Preferred Frame Theories 

Suppose we have a theory in which light and gravity propagate at fixed speeds, but in which 
c 7^ Cg. Such a theory will have two sets of "null cones," one determined by the propagation of 
gravity and one by the propagation of light, and thus two metrics.* It will also violate Lorentz 
invariance, in the sense that it will admit a preferred frame, the frame in which the speed 
of gravity and the speed of light arc both isotropic. Such violations are strongly constrained 
by experiment, but the strongest constraints are nongravitational [15]; as long as light and 
matter see the same metric, there is room for the gravitational metric to be different. 

The simplest, although certainly not the only, way to build such a theory is to start with 
a gravitational metric Qah and a unit timelike one-form field Ua, and to introduce a second 
metric Qab to describe the propagation of light and matter, 

9ab = 9ab - eUaUh, Q^^UaUb = 1, (2.1) 

where e is a fixed parameter. To keep notation simple, indices will be raised and lowered by 
Qab only; in particular, g""'^ = g°''^gcd9^'^ will not be the inverse of cjab- I choose coordinates 
such that gab reduces to the standard Minkowski metric rjab in fiat spacetime; in the preferred 
frame Ui — 0, we then have gab = diag(l — e, —1, —1, —1). The speed of light is thus 

c = vT^ (2.2) 

in units Cg — 1. If c > Cg, the absence of gravitational Cerenkov radiation places extremely 
strong limits on the difference [16]; I therefore take e to be positive. 

Let us now consider a theory with three dynamical sectors: gravity, the vector field u, and 
matter (including light): 

1. Gravity: I take the gravitational action to be the usual Einstein-Hilbert action for the 
metric gab, 

Igrav ^ ^ j (2-3) 

In standard general relativity, the coupling constant k depends on c; here, that dependence 
will be determined below from the Newtonian limit. 



*There are actually a number of independent ways in which the parameter c can enter; see [14]. I consider 
only the simplest. 
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2. Vector: The vector field Ua has four possible kinetic terms quadratic in derivatives [17]. I 
choose the simplest one, that of [12]: 

4 = ^ y d'^x^ [VaUbiV^u' - V'u'') + Xig'^'uaUh - 1)] . (2.4) 

With this choice the longitudinal mode becomes nonpropagating, and the propagating modes 
of u and g all travel at speed Cg] thus the "speed of gravity" is unique and well-defined [17]. 
The Lagrange multiplier A ensures that u remains a unit timelike vector. The coupling con- 
stant (3 is, for now, arbitrary; we shall see later that a particularly interesting case is that in 
which it is of order e. 



3. Matter: As noted above, nongravitational observations strongly suggest that all forms of 
matter see very nearly the same "speed of light." I therefore take a generic matter action to 
be of the form 

Imat ^ ImatbP,g] (2.5) 

with the standard minimal coupling to g. The definition of the stress-energy tensor in a 
bimetric theory is ambiguous; I choose the convention 

SImat d^xV^T''%ab ^\ j d^x^ {T'^'Sgab - 2eT^'ubSua) , (2.6) 

where the second equality comes from the fact that I^at depends on g and u only through 
the combination g. Note the choice of rather than y/—g in the integral. A different 

choice would introduce extra factors of c/cg. Once the final results are expressed in terms of 
physically measured quantities, however, this ambiguity will disappear. 

Two special cases will be important. The first is that of a point particle of mass m moving 
along a world line 7, 

/f f - /_ dx"' dx^\^^'^ 

ds^mj d'^x J dXS^{x - z{X)) \ 9ab-^-^J , (2.7) 

where 5 is a "densitized" delta function. An easy computation gives 

T^ = m / ds5\x - z{s)){l - e{u ■ vfY^'^v'^v^ with v" = — (2.8) 

where 5 is the ordinary ( "undensitized" ) delta function, i.e., J d'^x^J—g 5^{x) — 1. 

The second special case is light, which we can describe by the standard Maxwell action 
coupled to the metric g, 



I- EM 



^ j d^x^gr~9'\daA^ - d,A:){d,Aa - d^A,). (2.9) 



By standard geometric optics arguments (see, for example, [18]), it follows that light will 
follow the null geodesies of the metric g. To compute Shapiro time delay, we will thus need 
to understand the dynamics of this metric. 
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3 Field Equations and the Weak Field Approximation 

The field equations are easily obtained from eqns. (2.3-2.6): 

G"^ = {T^ + T^^) - pXu'^u'' (3.1) 

V„ (VV - VV) = -jeT^'ua + Xu" (3.2) 

g'^'uaUk = 1 (3.3) 

where 

= - A (^H\H'"' - ^g"-^H^dH'^^ with H"-^ = VV - VV. (3.4) 

Differentiating (3.1) and (3.2) and using the identities VaG"^ = and VoShH"^ = 0, we 
obtain a conservation law 

V„T«'' = eV„ [T^'u^u') - eT-«„V V. (3.5) 

It is not hard to show that this implies that VoT^ — 0, where V is the covariant deriva- 
tive compatible with the metric g. This means, in particular, that the field equations are 
compatible with the matter equations of motion in the geometry determined by g. 



3.1 Weak field approximation 

We next need the weak field approximation of these field equations. Let 

9ab = Vab + hab: = «a + Wa (3.6) 

where w is a constant timelike one- form such that r]°'^UaUb — 1, and where h and w are small. 
Define 

/ = da (^h^' - ^v'^'h^ , (3.7) 

where, as usual in the weak field approximation, we now raise and lower indices with the flat 
metric rj. It is then straightforward to check that to lowest order, the field equations become 

OKb - daXb - dbXa = -2«' (Tab " ^VabT^ + 2/3A (uaUb - ^Vab^ (3.8) 
nwb = dbidaw'') - —eTbcu" + Xub (3.9) 

2W„W" - U'^UaUb = 0, (3.10) 
where T"^ is now the matter metric alone {T^ vanishes at this order) . 
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The "matter" metric g can similarly be expanded: 

9ab = fjab + Kb ~ iVab - ^UaUh) + {Jlab - eUaWb - eUbWa). (S-H) 

Combining (3.8) and (3.9) and choosing a gauge 

Xa^Xa- euadbw'' = 0, (3.12) 

we find that 

Dhab = -2«;^ (Xab - ^VabT - 1^ {UaTb^Uc + UbTa^Uc)^ + 2(/3 - e)XUaUb - (3Xriab- (3.13) 

Note that if /3 = e <^ 1, this has nearly the same form as Kopeikin's modified field equations, 
eqns. (2.9) and (2.16) of [5]. 

To proceed further, we also need to understand the dynamics of A. Note first from (3.2) 
that A is small (of order w). Taking the divergence of (3.2), we then see that to the order at 
which we are working, 

u^\7^X = 0, (3.14) 

i.e., A is nonpropagating [17]. More precisely, choose a coordinate system u = (1,0,0,0), so 
that by (3.10), — ^hP^. Decompose into a transverse component {diW^'^ — 0) and 
a longitudinal component d^u;. Then (3.9) yields 

A = V^douj - ^hoo) + jeToo (3.15) 

with no further equation for cu. We thus can — and for the remainder of this paper, will — choose 
initial data for u such that A = at the required order. 

(One might worry that such data are unstable, since for a static source, (3.15) implies that 
u grows linearly in time. This instability is an artifact of our approximation, however, and 
disappears at higher orders [12]. Note that in the vacuum case, A = for the analog of the 
Schwarzschild metric [12].) 

3.2 Newtonian limit 

Let us next consider the Newtonian limit. Since matter moves along geodesies with respect 
to the metric g, this limit takes the usual form 

d^x* 1 ~ 1 ~ 

~ ^g'^djhoQ ^ --dihoQ, (3.16) 

where the last approximation requires that the velocity with respect to the preferred frame 
be much less than c, and therefore negligible in this limit. For the correct Newtonian limit, 
we thus require that /iqo — 2$, where $ is the Newtonian potential. By (2.8) and (3.13), on 
the other hand, 

V%o ~ 2k' (too -\t- ^Too^ ^^(^- (1 - t)-"'m5\x - z). (3.17) 



5 



We thus find that 



2e 



(3.18) 



4 Shapiro Time Delay 

We are at last in a position to compute the Shapiro time delay for this class of theories. 
We follow the basic procedure of [3]. In the absence of gravity, light follows a straight line. 
In the presence of a weak field, a light path can therefore be written as 



with small. (The suffix e denotes the point of emission.) Let 

dt dt 



(4.1) 



(4.2) 



be the tangent vector to the light's trajectory, with = 1 and = 0. As noted in section 2, 
this vector must be null with respect to the metric gab = Vab + hab- Thus to lowest order 



iVab - eUaUb)k''k^ = = 1 



e{uQ - k • u)^ f=i (1 - e) - |kp + 2ek • u. 



(4.3) 



where in the last line I have dropped terms of order |up. Not surprisingly, the speed of light 
is anisotropic: the speed in the k direction is 



Ck = (1 - e + 2ek • u) 



,-,^V2 



c + e- 



ku 



(4.4) 



At the next order. 



dv^ 

'^{r]ab - eUaUb)k''^ + habk^k"" = 0, 
dt 



implying, again up to terms of order |up, that 



d 



d 



^(k•y) = -e^(u•y) + ^/^„,rA;^ 



dt 

Now consider the square of (4.1): 

|2 . 



dt 



cl{t-t,y + 2k-y{t-te), 



or, using (4.6), 

t-te 



X — 



Ck 



k • y |x - Xe 



Ck 
|x - Xp 



4/ 

•^k Jte 



dt 



2 



+ e 



u-y(t) 



Ck 



+ A{t,t,), 



(4.5) 



(4.6) 



(4.7) 



(4.8) 
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where to this order the integral is along the unperturbed straight light path. For e = 0, 
this reduces to the standard Shapiro time delay formula (note that my signature conventions 
differ from [3]). For e ^ 0, it contains two new features, anisotropic "preferred frame" terms 
proportional to k ■ u and e-dependent terms that might distinguish the speed of light and the 
speed of gravity. The former are certainly of interest, and might be useful for constraining 
models of this sort, but they are not to the issue being addressed in this paper. Therefore, 
following Will [3] and Kopeikin [5], I will discard these terms, by choosing coordinates vP — 1, 
u = 0. 

By (2.8) and (3.13), the contribution to (4.8) of a point mass m and four- velocity v"" is 
given by 



Dihabk'^k') = -2k' [ T^k'^k' --T- -Toak'' 



(4.9) 



= -IGTrGm (^1 - ^ j j ds 5\x - z(s)) 



2 1 



— T v-k 

2/3 



where 1 have used (3.18) for /t^ and discarded terms of order |vp. The solution takes the 
standard Lienard- Wei chert form 



habk'^k'' = (1 + 7)G'm- 



l-(2 + C)^ 



X - z{sr) \ - v{sr) ■ (X - z{sr)) 



where sr is the retarded time,''^ 
and 



sr = t 



HSr 



(4.10) 
(4.11) 



2e2 

1 + 7-2(1-- 



2 P 



l-(l-6) 1- 



2(5) 



(4.12) 



I have chosen notation in such a way that this expression is identical to that of [3]; one can 
therefore read off the results directly from that paper (although one must keep careful track 
of factors of c) : 



A(i,ie) 



'j)Gm 



1-(1 + C)^^ ln(r-K-r^ 



with K = k k X (va X k), k = k/c. 



(4.13) 



There is one basic physical difference between this expression and the results of [3], though: 
while the expression in [3] had no apparent dependence on the speed of gravity, here such a 
dependence is explicit. 

^Note that since □ in (4.9) is the gravitational d'Alembertian, this expression depends on Cg, not c. 
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Two special cases are of interest. First, suppose e is small and (3 is of order unity. Then 
1 + 7 ^ 2(1 — e/2) ^ 2{c/cg) and 1 + C~1 — e = [cjcgf' . The first-order time delay 
thus acquires a prefactor Gm/cgC^, while the velocity-dependent prefactor is proportional to 
k-v/cg^. To allow more general possibilities, suppose instead that /? is of order e, say /5 = e/6, 
with e small. Then 

l + 7«2(l-(l-26)|)«2(c/c,)^-^'' 

1 + C ~ (1 - (1 - h)e) ^ {c/cgf^^-''^ (4.14) 



and thus 

A(t,te 

with 



2^771 



C2^ 



In(r-K-r)^ (4.15) 



2(l+i)) 1-26 >, 1 . ^ 

Cl = c 3 C2 = c''c/-^ (4.16) 

We see that by varying the coupling /3, one can make the velocity-dependent prefactor in 
(4.13) depend almost arbitrarily on c or Cg. The logarithmic term, on the other hand, depends 
only on the speed of light c. Since in ordinary general relativity this term gives the dominant 
contribution to the Fomalont-Kopeikin measurement, a bit more detail may be in order. The 
logarithm comes from an integral of a retarded potential, 

(4.17) 



z(s,^)|-^M^z£M' 



where I have restored the exphcit factor of Cg. The key observation — explained in [19] for the 
case of gravity and in [20] for the case of cicctromagnetism, and used by Will [3] — is that to 
first order in f , the denominator of this integral is independent of Cg. This does not, of course, 
mean that gravity does not propagate at the speed Cg. But to first order, the measurable 
effects of the propagation delay are canceled by velocity-dependent terms in the interaction, 
and observations are insensitive to the value of Cg in (4.17). More precisely, using (4.11), one 
easily finds that the denominator of (4.17) is 

|x - z(.^)| - ^M-(^-^M) ^ |x - z(.«) -{t- sn)^{t)\ + Oiv'/cg") 

= |x - z{t)\ + 0(^;7c/) = |x(t,) - z{tr) + (k - v)(t - i.)| + ^(^^Vc/) (4.18) 



where t is the "instantaneous" time and tr is the time of reception of the ray. Observe that 
the speed of gravity Cg has dropped out of the expression, which depends, at order v/cg, only 
on the unretarded "instantaneous" position of the source. Thus to this order, the integral 
(4.17) cannot depend on Cg} 

An explicit evaluation of (4.17) confirms this conclusion. The logarithmic term in (4.13) 
contains factors of c not because the retarded potential depends on c„ — this factor has canceled 



■I- This is the underlying reason that Faber's computation with Cg = oo [8] agrees with that of Will [3]. 
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out at order v/cg — but only because the integrand depends on |k — v|, and thus on |k|. By 
(4.3), |k| = c, and by (4.1) the c appearing here is directly the speed of Hght, that is, the 
distance traveled by the light ray divided by its travel time. As a straightforward check of this 
result, one can repeat the same computation for a rapidly moving massive particle traveling 
at a subluminal speed w near c. It is not hard to show that to first order in the only change 
in the argument of the logarithm in (4.13) is that v/cis replaced by v/w. 

5 Conclusion 

In ordinary general relativity, the speed of gravity and the speed of light are equal. This 
is not merely a coincidence; the two speeds express the same information about the metric 
structure of spacetime. Most of the controversy over the interpretation of the Fomalont- 
Kopeikin observation has implicitly assumed that there is some unique way to separate these 
speeds, that is, some unique deformation of general relativity to a c 7^ theory. We have 
now seen explicitly that this is not the case: we have constructed a class of theories that 
all reduce to general relativity as c — > Cg, but in which the Shapiro time delay depends on 
different combinations of c and Cg determined by the choice of an otherwise arbitrary coupling 
constant. 

In the process, we have reproduced one result of Will [3]: the vector K in (4.13) depends 
only on c, essentially because the c^-dependence of a retarded potential cancels to lowest order 
in w/c [19]. This is a generic feature of retarded potentials, and it would take a considerably 
more drastic shift away from general relativity to change this result. We have also, however, 
reproduced a result of Kopeikin [4,5]: if c 7^ Cg, conservation requires a modification of the 
stress-energy tensor, and thus a new Cg dependence in the equations of motion. Together, 
these results make the c- and c^-dependence of the Shapiro time delay complicated and highly 
model-dependent. The question of whether the time delay depends on the speed of gravity is 
thus unanswerable: without a good deal of further physical input into the choice of a c 7^ 
model, the dependence simply cannot be specified. 

For the models I have studied here, the best existing observations restrict 7 to I7 — 1| < 
5 X 10"^ [21]. By (4.14), this requires (1 - 26)e < 5 x 10"^. Note that this does not in 
itself restrict c/c^; if b is close enough to 1/2, e can be large. Within the present framework, 
though — and given the observed restrictions on 7, and thus the "closeness" to standard general 
relativity — it can be shown that the Fomalont- Kopeikin measurement is mainly sensitive to 
the logarithmic term in (4.13).^ As discussed above, to lowest order in v this term depends 
only on the speed of light c. One can thus argue that within the particular class of theories I 
have investigated, present measurements are not yet accurate enough to determine Cg. 

§This is unintuitive, since one expects a logarithm to be less sensitive than a power. But VLBI measures 
differences in arrival times, thus essentially differentiating a logarithm whose argument is small. The relevant 
computation in the present model is identical to that within general relativity, for which details can be found 
in section 4 of [7] ; the velocity-dependent part of the argument of the logarithm gives the last term in equation 
(12) of that paper, which is large because of the smallness of the angle 6 separating Jupiter and the quasar. 
Note that because Cg > c [16], C2 in (4.15) must be either greater than c (for b k, 1/2) or approximately equal 
to c (for e « 0), so the prefactor cannot become much larger than in standard general relativity. 
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On the other hand, this analysis has also made it clear that even this weaker claim depends 
on the class of "comparison theories" under consideration. It is not easy to see how to change 
the c-dcpcndence of the argument of the logarithm in (4.13); the absence of Cg to first order in 
V follows from very general properties of retarded potentials. But it is not at all clear that this 
term will always be the most important one. For example, it is evident from (4.12) that one 
can choose e and /? in such a way that 1 + C becomes very large, thus allowing a c^-dependent 
prefactor in (4.13) to dominate over the logarithmic term. For the simple models considered 
here, such a choice would make 7 unacccptably large. But by adding additional terms to the 
Lagrangian (2.4) for m, one can easily introduce new couplings that can allow 7 and 1 + ^ to 
vary more independently. 

The particular models I have considered here may be testable by other observations. In 
particular, they have preferred frame effects that might, for example, lead to a prediction 
of anomalous tides [22]. This is a generic feature of models with c ^ Cg. there is a unique 
frame in which light and gravity both propagate isotropically, and motion with respect to this 
frame is, in principle, observable. On the other hand, even in the general class of theories 
considered here, I have examined only a narrow set. The general kinetic term for u has four 
coupling constants, three fundamental propagation speeds, and dispersion relations that make 
these speeds energy-dependent [17], and it might be possible to tune the couplings to make 
preferred frame effects small while maintaining a nontrivial expression for Shapiro time delay. 
This is an interesting question, but it lies outside the scope of this paper. 
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